clear
x = -5:0.1:5;
y = -5:0.1:5;
[x,y] = meshgrid(x,y);
z = zeros(size(x));
[n,~] = size(x);
for i=1:n
    z(:,i) = mvnpdf([x(:,i) y(:,i)]);
end
mesh(x,y,z)
